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ABSTRACT:- A biorthogonal wavelet packets associated with nonuniform multiresoltion analysis 
(NUMRA) was introduced by Firdous shah. In this paper we generalize and define the biorthogonal 
nonuniform multiwavelet packets associated by nonuniform multiresolution analysis with multiplicity D 
(NUMRA-D). Further from the meaning of Fourier transform we study their characteristics. 
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I. INTRODUCTION 

In [9,10] Gabardo and Nashed considered a generalization of Mallat’s[16] theory of MRA based on 
spectral pairs where the translation set is of the form (0, r/N} + 2 Z, where N > 1 is an integer, 1 < r < 2N — 
1, r is an odd integer and r, N are relatively prime, which is called NUMRA. In [14] we provided the necessary 
and sufficient condition for the existence of NUMRAMW multiplicity D. G.Gripenberg and X. Wang gave the 
characterization theorem for dyadic orthonormal wavelet in L 2 (M) Wavelet packets constructed by coifman [15] 
it was gave rise to a large class of orthonormal bases of L 2 (M) into direct sum of its closed subspaces Wavelet 
packets. Behera[2,3] has constructed nonuniform wavelet packets associated with NUMRA. In [1] we 
constructed multiwavelet packets associated with nonuniform multiresolution analysis with multiplicity D 
(NUMRA-D). Firdous [8] introduce the notion of biorthogonal wavelet packets associated with nonuniform 
multiresolution analysis. Our result extends the biorthogonal wavelet packets of Firdous shah into NUMRA-D 
in has constructed nonuniform wavelet packets associated with NUMRA. 

II. BASIC DEFINITIONS AND NOTATION 

We will state some important preliminaries and notation in this section that we are need it in the recent 
paper. First we recall the definition of NUMRA — D as defined in [14], some of its properties and the associate 
multiwavelet packets [1-3, 8, 14] as follows: 

Definition 2.1. A nonuniform multiresolution analysis with multiplicity D for dilation 6 = 2Na and translation 
A, is a collection { V}j B i of closed subspaces of L 2 (R) satisfying the following axioms: 

(PI) Vj c V j+1 , for all j £ Z, 

(P2)/(-) £ Vj if and only if f (0 •) £ 1 j +1 , for all j £ Z, 

(P3) n ia Vj = {o}, 

(P 4) U ; 6Z Vj is dense in L 2 (W), and 

(P5) There exist functions r/j 1 , <p 2 , , (p D £ V 0 , called the scaling functions, such that the collection { <p d ( ■ 
— X) ■ X £ A, 1 < d < D} is a complete orthonormal basis for V 0 . 

In the axiom (P 5), the set of scaling functions <I> = { r/j 1 , c/j 2 , .... r/j ,J } is called multiscaling function of 
multiplicity D. When N > 1, the dilation factor of 6 ensures that PA c 2Z c A. 

In the theory of multresolution analysis, another sequence {Wj}j € i of closed subspaces of L l (M) is defined by 
Wj = Vj + i QVj , j £ Z and Q denotes the orthogonal complement of Vj in Vj +1 , for an NUMRA — D (Vy } 7 - 
with dilation factor 9. These subspaces hold the scaling property of (U}, eZ , and we have : 

L 2 (R) = © ;ez Wj = F 0 ®(© ; >o^) (1) 
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A set of functions {xp l ■ 1 < l < (9— 1)0} : = 'F in L l (M) is said to be a nonuniform multiwavelet 
associated with the NUMRA — D if the collection {?//(■ — A) : 1 < l < (0 — 1 )0, A £ A} forms an 

orthonormal basis for W 0 . We call 'F to be an NUMRA — D multiwavelet. In view of properties of 14} , the 
collection 

{■ 9l 12 ip 1 (9 1 ■ -A) : ; £ 1 , 1 < l < (0 - 1 )D, A £ A} 

forms an orthonormal basis for L 2 (M) if 'F is a nonuniform multi wavelet. 

Below we mention a result obtained in [14] that will be used in sequel. 

Theorem 2.2. Suppose {Fjyez is an NUMRA — D with dilation 9 and translation A. If there exist L functions 
xp k , 1 < k < L, in V 1 such that the family of functions 

V[ '■= {<p d {- - A), xp k f - A): 1 < d < D; 1 < k < L; A £ A} 

forms an orthonormal system for the generating subspace V 1 , then; L = (9 — 1 )D is a necessary and sufficient 
condition such that the above system is complete in V x . 



Proposition 2.3. Suppose {VJ- }y is an NUMRA — D with dilation 9 and translation A. Then the space V 1 
consists precisely of the functions / £ L 2 (R) whose Fourier transform can be written as 

D 



fm = Y J mf ' d co ^(o. 



a.e. / £ 



d = 1 



for locally L 2 functions mf ,d for 1 < d < D together with 



2ni$r 



( 2 ) 



(3) 



mf ,d {£) — + e n m'f (f), 

where m[' d and mf' 1 , d = 1, 2 are locally L 2 , 1/2-periodic functions. 

In the above proposition the space V 1 consist functions / £ L Z (R). Since (1 /9)f(x/9) £ V 0 , there exists a 

c f f f > 1 /\2 f ,D n-\ .• r • v -1 I f I 

sequence {a A = (aj ,a A ,...,a A '* N N 1 



)}i6A satisfying Sd=i Saga|^ 

l f Q = TY a2 ' d (pd{ - x ~ A) ’ 



< oo such that 



(4) 



d=l leA 

or, equivalently, by taking the Fourier transform of both sides of the previous equation, we obtained the result 
(2.1) and (2.2), with 



mf‘ d (f) = Y,x^a{' d e- 2 ^ 2 , (5) 

™{ 4 (f) = £ a{ 4 e~^ m , and m[ 4 (f) = £ a f £ +r _ e~^ m . (6) 

meZ rnETL 



We denote by O and *F column vectors in C D and in C L as ® = {( p^,...,(p D } and *F = {ip 1 ,..., ip L ], 
respectively. In particular, since <p d (x ) £ V 0 c V ] , from Proposition 2.3 there are locally L 2 functions mj’ , for 
1 < d,d < D, such that, for a.e. 

D D 



<, P d (9 f) = ^ m 0 

d'= 1 






(0 <?' 



® - I 

d'=l 

(p d ,d 



m. 



d,d 

01 



(0 + C 



2ni^r 



<2 (0)</(0 



(7) 



Taking = mj^ ' and = ^02 > in the matrix notation, this can be written as follows: 

0(00 = M 0 (O 0(0, for a.e. f E R, 



with M 0 (f) = (M 01 (f) + e ~ 2nl ^ r /N M 02 (/)) 

, where the matrices M 01 and M 02 defined by 

M 01 = ) and M 0 2 = [tnif ) , are usually called low-pass filters (or scaling matrix filters) 

' ’ \<fd,d <D ' ' 1 <d,d <D 

associated with the scaling family <f>. 

Similarly, in case of xp l ( x ) £ W 0 c F l7 where W 0 — span{ xp l (x — A): 1 < l < L; A £ A}, there are locally L 2 
functions mf ,d , for 1 < d < D and 1 < l < L such that, for a.e. £ R, 

D D 

V (00 = 2 (0^(0 = 2 (O + e~ 2ni ^ r ! N mf 2 d (f)) (8) 

d=l d=l 
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whose matrix notation is given as follows by considering m l {f = ,d and 



m 



W _ . 



12 



ra- 



il 



'PCef) = Mi (f) O(f), for a.e. f £ R, 

withMjff) = M n ( /) + e~ 2m ^ r /N M 12 (f), where the matrices M n and M 12 defined by 

Mu = «) Issusfs0 

and 

M 12 = ( m 12)i<;<L,i<d<D 

are usually called high-pass filters associated with the scaling family (T, <!>). 



Definition 2.4. The function of : 1 < l < L, 1 < d < D as defined as follows 

Sf(?) = lLi feR. (9) 

will be called the basic nonuniform multiwavelet packets corresponding to the NUMRA — D {Vj '■ j £ If of 
L Z (R) associated with the dilation 9. 

Note that (4.2) defines of for every non-negative integer 1 and every d such that 1 < d < D. 

Now Wg = (p d , d = 1 ,..,D is a multiscaling function of multiplicity D while of , 1 < l < 9 — 1, d = 1,- • • 
,D, are the basic multiwavelets associated with NUMRA — D. Now for any n £ N 0 = MU{0}, we define 
of:l < d < D, recursively as follows. Suppose that {of: 1 < d < D, p £ N 0 ] are defined already. Then 
define of +8p ,0 < q < 9 — 1, by 

Mq+ep W = 2Li ^eA 9 af q,h of ( 9x - A). (10) 

Note that (4.2) defines of for n > 0. Taking Fourier transform we get 

(tf +dp Xf) = SLi^'\f/0) of(f/9), 0 < q < 9 - 1, 1 < d < D. (11) 



Theorem 2.5. [1] Let [of :n> 0, 1 < d < D] be the basic nonuniform multiwavelet packets associated with 
the NUMRA - D {Vj}. Then 

(i) {of (■ —A ): ( \9y < n < (0) 7+1 — 1, A £ A, 1 < d < D }is an orthonormal basis of Wj , j > 0. 

(ii) {of (■ —A) : 0 < n < (0) ; — 1, A £ A , 1 < d < D] is an orthonormal basis of Vj , j > 0. 

(iii) {of (• —A) : n > 0, A £ A, 1 < d < D} is an orthonormal basis of L Z (R). 



Now as in (4), (7), (8), (9), (10) and (1 1) we can define: 

^m = i. D d ' = imi 4 \op d '(io 

V (flf)=Z2=i mf’ d (f)^(f) 

Sf (O = SLi ™ d ’ l ' h (9 v' 1 (f /8) ,a.e. / £ R 
tiq+ep W = 2Ll a f’ q ' h S p( 0X - 

(ti d q+ e v m = YLi^ d ' q ' h ^/9) 5*(f/0). 0 < q <9 - 1,1 < d < D 
And 



d,q,h 

m i H = 



/or q ^ 0; k 



'■ C\ + (u 



( 12 ) 

(13) 

(14) 

(15) 

(16) 



III. BIORTHOGONAL MULTIWAVELET PACKETS 
Definition 3 . 1 . A pair of function f(x), f(x) £ L z (E) are biorthogonal, if their translates satisfy 

</C),/(.-A)> = Scu- l£A. 

Where 5 0 2 is kronecker symbol. 

If cp d (A), <p d (x) £ L z (E) are a pair of biorthogonal scaling functions, then 

(<p d (.l<p d (.-X)) = 8 0 ' X , XeA,l< d<D. (17) 

Moreover we say that t pf(-), xpf(.—X)eL 2 {M.),l<r<9 — l,l<d<D are pair of biorthogonal 
nonuniform multiwavelet associated with a pair of biorthogonal scaling functions cp d (. ), tp d (. —A) if the family 
{ipf (. —A): A £ A, r = 1, ... ,9 — 1,1 < d < D] is a riesz of w 0 , and 
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{cp d f-X)) = 0, A E A,1 < r < 9 - 1,1 < d < D (18) 

{(p d {.),xpff-X)) = 0, Ae A,l<r <B - 1,1 < d < D (19) 

(xp? (. ), xp d (. -X)) = S rs S ojL , Ae A,1 <r,s <9 - 1,1 < d < D (20) 



Proposition 3.2. If xpf (xfxpf (x) E L 2 (R), 1 < r < 9 — 1 are a pair of biorthogonal nonuniform 
multiwavelet associated with a pair of biorthogonal scaling functions cp d (x), q> d (x), then L 2 (R) = 
©;a Wj = P 0 e(©;>o Wj) =@ je i ® e r Z\ W[ 

Where 

Wf = span [ip d ((9)1 . —A) : A E A,j E TL, 1 < r < 9 - 1, d = 1, ...,D] (21) 

Lemma 3.3. Let <p rf (x), <p d (x) be a pair of scaling functions. Then (p d (x), <p d (x) are biorthogonal scaling if 
and only if £d=i £ AeA <p d - A) q> d {f - X) = 1 a. e % E R. 

Lemma 3.4. Asume that o> d , 3 d E if (R), 1 < q < 0 — 1, 1 < d < D are a pair of biorthogonal nonuniform 
multiwavelet associated with a pair of biorthogonal scaling functions nig (x) and Sg (x). Then, 

Z®=o ZLi md ' p ' h (f + 2na ) ™- i,q,h + 2ncj ) = <W Q < p,q < 9 — 1,1 < d < D (22) 

Proof. By (11), (14), (17), (21) and Lemma 3.3, we have 

Sp, q = ^ «Jp (? + 2ttX) 3 d (f + 2 nX) 

AgA 

D 

= X 2 md ' P ' h ( 9 ~ 1 ^ + 2?r '0) s o (0 -1 (f + 27r2)) 

16A/i=1 

x + 2 tt/1)) + 2ttX)) 

0-1 D 

= ^ m d,p ' h (9~ 1 (f + 27 rcr)) rh d ' q ' h (^9~ 1 (f +27 rcr)) 

(T — 0 /l = l 

X Sg (0 _1 (f + 27 T<t) + 27T/1}) 3q (0 _1 (f + 27T(j) + 27T/1) 

9-1 D 

= ^ m d,p,h {9~ 1 {f + 27rcr)) m c! '' ? ' ,l (0 _1 (f + 27 to-}) . ■ 

<7=0 /l=l 



Theorem 3.5. Suppose { u> d (x)\n > 0, 1 < d < D }and { 3 d (x):n > 0, 1 < d < D} are nonuniform 
multiwavelet Packets of multipilicity D with respect to a pair of biorthogonal scaling functions o)q (x) and 
3q (x). respectively. Then, for n > 0, we have 

{co d {.),3 d {.-X)) = 8 0A , A E A,1 < d < D (23) 

Proof. The proof it will be by induction, for n = 0,1, ... ,9 — 1, (15) is true from (11), (14). Assume that it 
holds when n < l, l > 0. For n — l. Order n — q + 9p , where p>0,0<q<9 — l&p<n. We have 
(o> d (. ), 3 d (. -A)) = S 0A o ZS =1 Z AeA 3 d (f - 2) 3 d (f-A) = 1,(61. 

And 

(wn (■), 3 d (.-/!)> = (3 d (. ), S d (.-A)) 



= f 8) d +gp (0 3 d +gp (f) e 2nU f df 

Jr 



= [ y° m d ' q ’ h (f/ 9 ) 3 d (f/ 9 ) y D m d ' qdl {f/ 9 )(o d {f/ 9 ) e 2ni ^ df 
Je ^ 1/1 = 1 

= yy D C md,q ' h (I + 2 ^) (f + 27 ta) m d ' q ' h (f/ 9 ) 5 *(f/ 0 ) e 2ni df 

,/l=1 ^([0,2 tt]+A) / \& ) 

( D \ / D \ 

I 1 (I + H 4 + 2,r;i) ) x (Z ‘ (I) 

AeAft=l / \h=i 



] e 2iri ^ df 



= f y m d,q ' h m d ’ q ' h (f/ 9 ') e 2ni ^ df 

^ro.277-i “ 



J 6 [0,27t] “ 
9-1 D 



= J ^ ^ m d ’ q ’ h + 27 ro-^ m d ’ q ’ h (— + 27 TO-) e 2m;i ^ d( 

l 0 ' 2 ''] ff=0 h= i 
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e 2nU t df = S 0il . • 

J[0,2n] 

Theorem 3.6. Suppose {a> d (xy.n > 0, 1 < d < D] and {co d (xy.n > 0,1 < d < D] are Nonuniform 
multiwavelet Packets of multipilicity D with respect to a pair of biorthogonal scaling functions 0 )q (x) and 
So (x), respectively. Then, for n > 0, we have 
<< +flp (■ X tf 2+Bp (■ -A)> = Sru S qi , q2 , A 6 A, 0 < q v q 2 < 9 - 1,1 < d < D (24) 

Proof. Similar steps in above proof. 

Theorem 3.7. If {a> d (x):n> 0,d — 1, ...,D} and {5) d (x):n> 0, d — l,...,D] are basic nonuniform 
multiwavelet packets with respect to a pair of biorthogonal multiscaling functions u)q(x) and 6) d (x), 
respectively. Then, for m, n > 0, we have 

«(.),S n d (.-A)> = S m ,„ S QJL (25) 

Proof. For m — n, the result (25) follows by Theorem 3.5. When m =4 n, and 0 < m, n < 9 — 1 0 < f, the result 
(25) can be get from Theorem 3.6. Assuming that m 4= n, and at least one of { m , n } does not lies in 0,1, ... ,6 — 
1, then we can rewrite m,n as m = q ± + 9p 1 ,n = q 2 + 9p 2 , where p 1 p 2 > 0, 0 < q lt q 2 < 9 — 1. 

Case 1. If Pi = p 2 , then q 1 =4 q 2 . Therefore, Eq. (25) follows by (11), (15), Lemma and (22), i.e.. 



« (. ), <4 (. -X)) = << +epi (. ), < +flp2 (. —X)) 



= <<+e P ,C),^ 2+flp2 (. —X)) 

= f S U,(')l^W)e 2 -^ df 

E 

= j m d,qi,h (9~ 1 O a&Cfl -1 O S d 2 (9- 1 O m d 'K' h ( 0" 1 0 e 2 * u f df 

E 

= ^ J m d ' qi ’ h (9- 1 0 f) m d ' q i’ h ( 0 _1 f) e 2irU ^ df 

^eAe([o,2jr]+l) 

= J (f + A)) S^Ce- 1 (f + 4)) ) m d ' qi,h ( 0 _1 f) fn d - <?2 ' ,l ( 0 _1 O e 2ni2 ^d^ 

' leA / 

0-1 

= J ^ m d,qi ’ h ( 0- 1 (f + 27TO-)) m d ’ q 2' ,! ( d” 1 (f + 2 tto-)) e 2 11 f 1 * df 
= J S qi , q2 e 2 ^ d(= S 0A = 0 

[0,2jt] 



S[0,2jt] Vie A 



[0,2tt] 



Case 2.U p p 2 , order p 3 = q 3 +0p 3 , p 2 = q 4 +0 p 4 , where > 0 and 0<q 3 , q 4 <0 — l.lf 

p 3 — p 4, then <74 4= ^4. Similar to Case 1, (25) can be established. When p 3 ^ p 4, we order p 3 — q 5 + 
0 p 5 . p 4 = + 0 p 6, where 5 /> 6 >0 and 0 < q 5 ,q 6 < 0 — 1. Thus, after taking finite steps (denoted by 

/f), we obtain 0 < p k ,p 2k ^0—1 and 0 < p k q 2k < 0 — l.\i p k — p 2k , then ^ =4 Similar to the 
Case 1, (25) follows. If p k =4 q k 
, then 

«, C),2w (■)> = 0 - ->O = 0 



Furthermore, we obtain 



X u n (. -a )) = (S OT (. ), ^ (. -y? )) 



^ d d 

V 7 + & P j (') ’ ^q2+&P2^-' 
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= I ^ l+ffpi {.)tf„ 2+0p2 {.-A) e 2na f d? 

R 



200) □ 0300-2000400-2000,02,A0-1000,04,A0-1002000000 



= ; R (n k r =i * P »w- k n (n %i n)e 2 ** # 



2k 



A£ W(\0,2ir}+A) ' r =l 2 



V=J 



l {V\^r.h {0 -r n \ /" Y[ ri**r + l*(0-r f) \ £ (f 



gk\0.2rr ) \ r =l 



\ r=l 



k£A 



+ Af)%„ 2 A0~\f + A))e 2 ”* d? 



= 0 m 
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